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I illustrate the role of infrared renormalons in computing inclusive B-decay spectra. 
I explain the relation between the leading ambiguity in the definition of Sudakov 
form factor ~ exp(JVA/M) and that of the pole mass, and show how these ambi- 
guities cancel out between the perturbative and non-perturbative components of 
the b-quark distribution in the meson. 

1. Introduction 

B-decay physics is gradually turning into a field of precision phenomenology. 
Inclusive decay measurements provide some of the most robust tests of the 
standard model. Classical examples are the rate of B — ► X s ^ decays [1] 
and constraints on the unitarity triangle through the measurement of V u b 
from charmless semileptonic decays [2]. 

The advantage of inclusive measurements over exclusive ones is that the 
corresponding theoretical predictions are, to large extent, free of hadronic 
uncertainties. QCD corrections to total decay rates are dominated by short 
distance scales, of order of the heavy-quark mass m, and are therefore 
primarily perturbative. Confinement effects appear as power corrections 
in A/m. Moreover, the Operator Product Expansion (OPE) allows one 
to estimate these power corrections by relating them to specific matrix 
elements of local operators between B-meson states, which are defined in 
the infinite-mass limit in the framework of the heavy-quark effective theory 
(HQET) [3]. These matrix elements can either be computed on the lattice 
or extracted from experimental data. 

In reality, however, experiments cannot perform completely inclusive 
measurements. Precise measurements are restricted to certain kinematic 
regions where the background is sufficiently low. The experimentally ac- 
cessible region in B — ► X s ^ is where the photon energy E 1 in the B rest 
frame is close to its maximal possible value, M/2, (M is the B meson mass), 
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or, equivalcntly, x = 2E~ f /M is near 1, which is the endpoint. Similarly, the 
accessible region in the CKM-suppressed B — ► X u lv decay is where the 
lepton energy fraction is near maximal, or where the invariant mass of the 
hadronic system is small. Out of this region this decay mode is completely 
overshadowed by the decay into charm. 

As a consequence, precision phenomenology must rely on detailed the- 
oretical understanding of the spectrum [4] . Of particular importance is the 
spectrum near the endpoint. It turn out, however, that the endpoint region 
is theoretically much harder to access as both the perturbative expansion 
and the OPE break down there. In the large- X region gluon emission is re- 
stricted to be soft or collinear to the light-quark jet. While the associated 
singularities cancel with virtual corrections (decay spectra being infrared 
and collinear safe) large Sudakov logarithms of (1 — x) appear in the expan- 
sion, which must therefore be resummed. Moreover, the OPE breaks down 
since the hierarchy between operators scaling with different powers of the 
mass is lost when (1 — x)M become as small as the QCD scale. Physically 
this reflects the fact that the spectrum in the endpoint region is driven by 
the dynamics of the light degrees-of-freedom in the meson. 

The lightcone-momentum distribution a of the b-quark in the B-meson 
has a particularly important role in the endpoint region [5-13]. It has been 
shown that up to subleading corrections 0(A/m) the physical spectrum 
can be obtained as the convolution between a perturbatively-calculable 
coefficient function and the QDF, where the latter essentially determines 
the shape for x — ► 1 . The key point is that the QDF is a property of the B 
meson, not of the particular decay mode considered, so it can be measured 
in one decay and used in another. Moreover, a systematic analysis of the 
QDF in the HQET highlights the significance of a few specific parameters 
which constitute the first few moments of this function: most importantly 
A = M — m, the difference between the meson mass and the quark pole 
mass, and then Ai corresponding to the kinetic energy of the b quark in 
the meson. 

Nevertheless, the dependence of theoretical predictions for the spectra 
on the QDF is still a major source of uncertainty. Apart from identify- 
ing its first few moments, very little is known about this function, so the 
phenomenology of decay spectra in the immediate vicinity of the endpoint 



a We shall define this function in full QCD, and call it Quark Distribution Function 
(QDF). This should be distinguished from the common practice to define it directly in 
the HQET, where the name "Shape Function" is often used. 
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(x — ► 1) remains, to large extent, model dependent. On the other hand, 
successful precision phenomenology can well be expected for more moder- 
ate (yet large) x values, corresponding to the region where the distribution 
peaks. Here the main obstacle has been in combining [10, 14] perturba- 
tive Sudakov effects with the HQET-based non-perturbative treatment dis- 
cussed above. 

It has recently been shown [15] that the resolution of this problem is 
firmly connected with infrared renormalons (for general review of renor- 
malons see [16]). Since the formulation of the HQET as well as the per- 
turbative calculation of decay spectra rely on the concept of an on-shcll 
heavy quark, both ingredients suffer from renormalon ambiguities. These 
ambiguities cancel out, of course, in the physical spectra. It is therefore 
useful to trace b the precise cancellation of ambiguities: the use of the 
HQET brings about dependence on the quark pole mass, which has a linear 
renormalon ambiguity [17-20]. This ambiguity cancels against the leading 
renormalon ambiguity in the Sudakov exponent [15]. In order to achieve 
power-like separation between perturbative and non-perturbative contribu- 
tions to decay spectra, one must therefore compute the Sudakov exponent 
as an asymptotic expansion, thus replacing the standard Sudakov resum- 
mation with fixed logarithmic accuracy by Dressed Gluon Exponentiation 
(DGE) [15,21-26]. 

In what follows we illustrate the role of renormalons in the QCD descrip- 
tion of decay spectra. We begin by briefly reviewing the HQET analysis for 
the QDF where we identify dependence on the quark pole mass. We recall 
that the pole mass suffers from an infrared renormalon ambiguity and show 
how this affects the QDF [15]. We then consider inclusive B-meson decays 
within perturbation theory, review the relevant results on large- a; factor- 
ization and Sudakov rcsummation [13], and then show that renormalon 
ambiguities appear in the Sudakov exponent [15], which, we emphasize, is 
a general phenomenon rather than a peculiarity of B decays. Finally, we 
combine the perturbative and non-perturbative ingredients recovering an 
unambiguous answer for the QDF in the meson and consequently for decay 
spectra. We conclude by shortly discussing the implications for precision 
phenomenology in inclusive decays. 



b This can be understood in analogy with factorization scale dependence, the main dif- 
ference being that here the interest is in power terms. 
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2. Heavy-quark effective theory and the QDF 

We define the QDF f(z;fi) as the Fourier transform of the forward hadronic 
matrix element of two heavy-quarks fields on the lightconc (y 2 = 0) : 

1 f°° + - 

f{z;n) = — J_Jy- C ^v- {B (p B )\^{Q) 1 +^(y)\B{ PB ))^ (1) 

where a path-ordered exponential between the fields is understood, is the 
B-meson momentum (p B — M 2 ), z is the fraction of the "+" momentum 
component carried by the b-quark field and fi is the renormalization scale 
of the operator. Decay spectra can be computed as a convolution between 
a perturbatively calculable coefficient function and f(z;fi). Let us first 
analyze f(z;n) non-perturbatively — we denote it f Nr (z) — suppressing 
any perturbative corrections. These will be recovered later on. 

Since the 6-quark mass is large, the heavy quark is not far from its mass 
shell. This observation is the basis of the HQET. The momentum of the 
heavy quark is p = mv + k where v is the hadron four velocity, v = p B /M, 
and k is a residual momentum, |fc| <C to. The effective field is defined by 
scaling out the dependence on the quark mass: h v (x) = e lmv ' x ^(1+^) ^(x). 
It then follows from the definition (1) that in the heavy-quark limit [5-15] 

J\zf^{z)^ tv - V(K - {1 - Z)M) = ^ (B(Mv)\h v (0)h v (y-)\B(Mv)) 

= 1 + ^(-iv ■ yf + ^(-iv ■ yf + • • • = H-iv ■ y), (2) 

where we inverted the Fourier transform and defined A = M — to; in the 
second line we expanded the lightcone operator in the HQET in terms of 
local operators, where e.g. 

-3/ 2 - A! = ^{B(Mv)\h v (0){ Sia ,-VpV v )iDi i iD v h v (0)\B(Mv)). (3) 

The HQET matrix elements /„ ~ 0(A n ) do not depend on the definition 
of the mass. On the other hand the vanishing of the linear term in (—iv ■ y), 
fi = 0, (and the absence of additional, mass dependent terms in front of 
higher powers of (—ivy)) in the second line of Eq. (2) are due to the HQET 
equation of motion for the heavy quark. Thus Eq. (2) relies on using the 
pole mass to define the HQET C . 



c The use of the pole mass in the field redefinition can be avoided if a residual mass term 
is introduced. This, however, does not change any of the conclusions [15]. 
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Expanding the exponential on the l.h.s of Eq. (2) we obtain: 

fn= [ /nf(*) (A-(l-z)M)", (4) 
Jo 

so the moments are fixed by the local matrix elements in the HQET. Mcllin 
moments are defined by 

K*= I dzz N -\f NP (z). (5) 
Jo 

For the first few Mcllin moments we have Ff p = f and 

pNP m . kiNP 

_/m\ 2 _lAi _fm\ 3 _X 1 m f 3 

2 ~M' 3 ~\m) 3M 2 ' 4 ~vmJ M3 + M3- W 
It is apparent that all the moments depend on the quark pole mass. They 
satisfy MdF^/dm = (N - 1) F^_ v At large N they are given by [15] 

= e-^ N -^ M ^{{N - 1)/M) + 0(1/AO, (7) 

where the exponential factor depends on the pole mass through A while 
T((N — 1)/M), defined in (2), is entirely quark-mass independent. Note 
that large N corresponds to asymptotically large lightcone separations. 



3. IR renormalon ambiguity in the pole mass 

The result of the previous section indicates inherent dependence of the 
non-perturbative component of the QDF in the heavy-quark limit on A = 
M — m, or, equivalently, on the pole mass m. We recall that the pole 
mass is defined in perturbation theory by requiring that the inverse quark 
propagator — mjg — S(p, m.Ms) vanishes at p 2 — m 2 . At any given order in 
a s one can solve the resulting equation obtaining a unique relation between 
the pole mass and the MS mass (or any other renormalized short-distance 
mass). However, when considered to power accuracy this definition remains 
ambiguous [17-20]. The on-shell condition brings about linear sensitivity to 
long-distance scales. In the perturbative expansion (in schemes such as MS) 
this sensitivity translates into non-alternating factorial divergence making 
the sum of the series ambiguous — the well known infrared renormalon. 

Specifically, with a single dressed gluon — thus to leading order in the 
largc-/? limit — the relation between the pole mass and ttjms is given by 
the following Borel sum [18]: 



Toms A) Jo \ m - 



1 

'•MS 



3eWl - u)T(u)T(l - 2u) 3 „ , , 

2r(3- u) Tu +R ^ u) 

(8) 
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where (u) is free of singularities and /3q = y^C/i — g-^Yf ■ The singularity 
of the integrand at u = \ translates into an 0(A/m) ambiguity, which 
directly affects the QDF moments F^ F ■ 

As one would expect, the pole-mass renormalon ambiguity cancels out 
whenever the pole mass is used to compute an observable quantity. A well- 
known example [19] is the calculation of the total semileptonic decay rate, 
which explicitly depends on the fifth power of the mass. Another example 
is the total energy of quarkonia [27]. Here we review this cancellation for 
the QDF in the meson and consequently for decay spectra [15]. 



4. Large-a; factorization in B decay 

Let us now consider B-decay spectra in perturbation theory. By taking the 
initial state to be an on-shell b quark we neglect non-perturbative effects as- 
sociated with the meson structure. The decay rate is infrared and collinear 
safe, so the partonic calculation yields finite perturbative expansion when 
expressed in terms of the renormalized coupling and mass. It should be 
noted, however, that this finitcness is owing to cancellation of logarithmic 
singularities between real and virtual corrections. As usual, the singularity 
leaves a trace in the form of Sudakov logarithms of (1 — x). Since such 
logarithms appear at any order in the perturbative expansion, they must 
be resummed. 

The resummation of Sudakov logarithms takes the form of exponenti- 
ation in Mellin space. This is a consequence of the factorization property 
of QCD matrix elements in the soft and collinear limits together with the 
factorization of phase space [13,29,30]. Up to (D(l/N) corrections the per- 
turbative expansion of inclusive decay spectra can be written in Mellin 
space as a product of three functions [13]: a soft function depending on 
m/N, a jet function depending on m 2 /N and a hard function depending 
on m — see Fig. 1. 

Furthermore, the resummation can be formulated as DGE, rcsumming 
running-coupling effects in the Sudakov exponent. The result is most conve- 
niently expressed as a Borel sum. Explicitly, for d B — ► X s j we have [15]: 

f 1 1 r/T PT 

M7= dx———x N - 1 = C N (m)J N (m ]fJl )S N (m;^ + 0(l/N), (9) 
Jo 1 tot dx 



d The corresponding formula for the semileptonic decay appears in [15]. 




where 



(C F f°° du , . /A 2 
bN{m;fJ,) = exp<^ — / — T(u) — 

B s (u)T(-2u) (N 2u - 



(10) 



— ) B A (u)\nN 



J N {m;n) = exp ^ — / — T(u) — 
I Po Jo u \ m 

- fl^(u)r(-v) (JV" - 1) - I ^ ) B^lniV 



(11) 



Here S^m;^) and J/v(m;/z) are the soft and jet functions, respectively . 
These functions were both normalized to unity — the exponents vanish at 
N = 1 — so they acquire dependence on the hard scale. Bs{u), Bj(u) and 
Ba(u) are Borel representations of anomalous dimensions of the soft, jet 
and cusp functions, respectively. In the large-/?o limit 



B s (u) = e c " (1 - u) +O(l//3 ), 

B ^-r {— + — 



o(VA) 



(12) 
(13) 



where c = 5/3 in MS. Beyond this limit the anomalous dimensions are 
known only as an expansion in u (through NNLO). Terms that are sub- 
leading in l/po, which appear first at 0(u r ) : are not small in QCD. The 
advantage of the large- /3q limit, where an analytic function is known, is that 
it allows one to verify the exact cancellation of renormalon ambiguities. 



e T(u) depends [15] on the approximation used for the /3 function; for one-loop running 
coupling T(u) = 1. 



Wilson line 

y- 




On shell b quark 



Figure 2. Process-independent calculation of the QDF in an on-shcll heavy quark in 
the large-/3o limit: the gluon is dressed by any number of fermion-loop insertions and 
then Nf > — 6/3o- In the A + = axial gauge only this diagram contributes. 

The soft function is the large- TV limit of the lightcone momentum dis- 
tribution of a b-quark field in an on-shcll b quark. In can be computed in 
a process- independent manner from the QDF definition (1), replacing the 
external meson states \B(Pb)) by on-shell quark states \b(p)). For example, 
the large-/?o Unfit result of Eqs. (10) and (12) can be obtained from the 
diagram of Fig. 2. It is related by crossing to the perturbative heavy-quark 
fragmentation function analyzed in [25]. 

The jet function describes the radiation associated with an unresolved 
jet of invariant mass m 2 /N. It is a universal object appearing in many ob- 
servables including deep inelastic structure functions [22-24], single-particle 
inclusive cross sections [22,25] and event-shape distribution [21,26,28]. 

In both the soft and jet functions there are renormalon ambiguities 
owing to the singularities of T(—2u) and r(— u), respectively. They appear 
as a result of integrating over the longitudinal momentum fraction near 
the endpoint and reflect the sensitivity of the exponent to large-distance 
scales through the running of the coupling. The ambiguity indicates the 
presence of non-perturbative power corrections at the exponent for each of 
the functions: powers of NA/m in the soft function and powers of NA 2 /m 2 
in the jet function. 

5. Cancellation of renormalon ambiguities in the exponent 

When considered to power accuracy, the perturbative soft function of 
Eq. (10) becomes ambiguous. This is not surprising since its definition in- 
volves the on-shell quark state \b(p)}. Its non-perturbative analog, the QDF 
in the meson defined in Eq. (1), should be well defined. Yet, at large N 
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these two functions differ just by (an infinite set of) non-perturbative power 
corrections on the scale M/N: 

F N {M;^ = S N {m;fj,)F^, (14) 

where f F^ p is given by (7). The corresponding non-perturbative large- a; 
factorization in B-meson decay is: 

f 1 1 dT 

M N = dx——x N - 1 =C N (m)J N (m; (J ,)F N (M;^ + 0(l/N), (15) 

*/ Q tot 

where the sole difference from the perturbative formula (9) is the replace- 
ment of the QDF in the quark, Spf(m; fj,), by that in the meson, F/v(M; y). 

Since the QDF F/v(M;/x) directly enters the measurable moments Mjv, 
it must be well defined. This will be the case only if renormalon ambiguities 
cancel in (14) between the perturbative and non-perturbative components. 
Indeed, such cancellation is expected because both Sjv(m;/z) and Fj^ p in- 
volve the concept of an on-shell heavy quark, while Fn(M;^j,) does not. 
Putting together (7) and (10) we obtain [15]: 

B s (u)T(-2u) (N 2u - 1) + (^pj B A {u) \nN J + 0(1/N), (16) 

which can be explicitly verified to be free of u = \ ambiguities in the 
large-/?o limit by substituting A = M — m and using Eqs. (8) and (12). 



6. Prospects for precision phenomenology 

The results of Sees. 4 and 5 have direct implications for the calculation of 
inclusive decay spectra: they open up the way for consistent power-like sepa- 
ration between perturbative and non-perturbative contributions depending 
on NA/M. QCD predictions for decay spectra in the peak region require 
Sudakov resummation as well as non-perturbative corrections depending 
on the meson structure. However, conventional Sudakov resummation with 
fixed logarithmic accuracy does not deal with the problem of separation 
between perturbative and non-perturbative contributions (the immediate 
price is Landau singularities). The perturbative coefficients get significant 
contributions from small momentum scales, contributions that increase with 



f We systematically neglect 0(1/N), or, equivalently, C(A/M) effects. Eq. (14) does not 
hold for small N. 
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increasing logarithmic accuracy [22]. DGE addressed the source of this 
problem. Using Borel summation — with a principal-value (PV) prescrip- 
tion for example — it systematically separates non-perturbative power-like 
terms in the Sudakov exponent from perturbative contributions. This pro- 
cedure uniquely defines the non-perturbative power terms — this is pre- 
cisely the meaning of Eq. (16): consider for simplicity the hypothetical 
case where T ~ 1 so f NP (z) ~ 5(z — m/M). The shape of the QDF is 
then determined by the perturbative Sudakov form factor; it is just shifted 
non-perturbatively by A/M toward smaller z values (in the general case T 
leads to some smearing). Precise control of this shift is, of course, crucial; 
an ambiguity of order A in A would be a catastrophe. However, based on 
Eq. (16) A is uniquely fixed: if the principal value of the Borel sum is used 
to define the perturbative Sudakov exponent, the same prescription must 
be used to relate the pole mass to any (well measured) short-distance mass 
when computing A, so A = M — m PV . 

It should be emphasized that quantitative control of power-like con- 
tributions by means of Borel summation requires more information than 
available either from fixed-order calculations of the anomalous dimensions 
or from the large-/3o limit alone. For example, the value of Bg(u) near u = \ 
becomes relevant. While challenging, this question can still be addressed 
within perturbation theory. 
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